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It has been recently pointed out that coupled dark matter-dark energy systems suffer from 
non-adiabatic instabilities at early times and large scales. We show how coupled models free 
from non-adiabatic instabilities can be identified as a function of a generic coupling Q and 
of the dark energy equation of state w. In our analysis, we do not refer to any particular 
cosmic field. We also confront a viable class of model in which the interaction is directly 
proportional to the dark energy density to recent cosmological data. In that framework, we 
show the correlations between the dark coupling and several cosmological parameters allowing 
to e.g. larger neutrino mass than in uncoupled models. 



1 Introduction 

Interactions between dark matter and dark energy are still allowed by observational data today. 
At the level of the background evolution equations, one can introduce a coupling between these 
two sectors in the following way: 



where PdmiPde) denotes the dark matter (dark energy) energy density, the dot indicates derivative 
with respect to conformal time dr = dt/a, TC = d/a and w = Pde/Pde is the dark-energy equation 
of state (P denotes the pressure). We work with the Friedman- Robertson- Walker (FRW) metric, 
assuming a flat universe and pressureless dark matter Wdm = Pdm/Pdm = 0. 

Q encodes the dark coupling and drives the energy exchange between dark matter and dark 
energy. For e.g. Q < the energy flows from dark matter to dark energy. It also changes the 
way that dark matter and dark energy redshift acting as an extra contribution to their effective 
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equation of state. In particular, for e.g. Q < 0, dark matter redshifts faster so that there is 
more dark matter in the past compared to uncoupled scenarios assuming that the dark matter 
density today is the same in the two models. Matter-radiation equality happens earlier and the 
growth of dark matter clustering is enhanced. This is one of the features which enable us to 
constraint the model with available cosmological data (see section [3|) . 

In order to deduce the evolution of density and velocity perturbations in coupled models, we 
need an expression of the energy transfer in terms of the stress-energy tensor. We follow 
in parameterizing the interaction as: 

Vm^Sh^ = Qui'-'^/a, (3) 
V.Tge), = -Q4'^™V«, (4) 
with TH n and T^, n the energy momentum tensors for the dark matter and dark energy 
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components, respectively. The dark matter four velocity ulf"^^ is defined in the synchronous 
gauge, in terms of the fluid proper velocity w^^^^, as uif'^'^ = a(— 1, v^^^^), where fi = 0..3 and 

1 = 1..3. This choice of parameterization guaranties the conservation of the total stress energy 
tensor of the system while it avoids momentum transfer in the rest frame of dark matter in 
which case one can work in the synchronous gauge comoving with dark matter (i.e v^,^^) = 0). 

We provide a criteria associated to the dubbed the doom factor (see section \2A\i to identify 
the stability region of coupled models satisfying to Eq. ([3]) and The doom factor is a function 
of the model parameters such as Q and w but it is defined independently of the explicit form of 
the coupling Q. Notice that the dark coupling terms which appears in the non-adiabatic dark 
energy pressure perturbations were first pointed out as a source for early time instabilities at 
large scales in Ref.^. 

Our results will then be illustrated within a successful class of models, in which Q is pro- 
portional to the dark energy density. Present data will be shown to allow for a sizeable in- 
teraction strength and to imply weaker cosmological limits on neutrino masses with respect to 
non- interacting scenarios. 

2 Origin of non-adiabatic instabilities 

Non-adiabatic instabilities arises at linear order in perturbations. Using the publicly available 
CAMB CO we preferred to work in the synchronous gauge comoving with dark matter. 

2.1 Doom factor 

Given our gauge choice, it is necessary to work out the expression for dark energy pressure 
perturbations SPde in the rest frame for dark matter. It can be shown that in the presence of a 
dark coupling (see e.g. Ref. ^ and also Re^5) it is given by: 

^^'^'^ 5Le + 3(cLe-cLe)(l + ^«)(l + d)^^, (5) 
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where 5pde denotes the dark energy energy density perturbation and 5de = ^ Pde I Pde-, ^de 
is the divergence of the fluid proper velocity v^i^de)^ ^Ide '^^ propagation speed of pressure 
fluctuations in the rest frame of dark energy and c^^^ = Pde/ Pde is the so called "adiabatic 
sound speed". In the following, we work with constant equation of state w in which case 
^ade ~ ^ ^^'^ assume that our universe is in accelerating expansion today which implies that 
w < —1/3. Moreover we restrict our analysis to the case c^^g > and c^de ~ ^ ^i^^ assumed 
for numerical computation. 



In Eq. ([5|) d refers to the doom factor which we have defined as: 
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We dub it so as it is precisely this extra factor, proportional to the dark coupling Q, which may 
induce non-adiabatic instabilities in the evolution of dark energy perturbations. Its sign will be 
determinant, as we show in the next section. 



2.2 Growth equation 



A cartoon equation of the growth equation governing the evolution of energy density linear 
perturbation for any species i,j is given by: 
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where 5, = bpij pi and ' = d/da. The evolution of a perturbation depends on the relative 
weight of the three terms present in the equation and on their signs: 

1. For positive A, the A and B terms taken by themselves would induce a rapid growth of the 
perturbation, which may be damped or antidamped (reinforced) depending on whether B 
is negative or positive, respectively §. In particular, for A and B both positive, the solution 
may enter in an exponentially growing, unstable, regime. 

2. For negative A, in contrast, the A and B terms taken alone describe a harmonic oscillator, 
with oscillations damped (antidamped) if B is negative (positive). In the A, S < regime, 
the third term may plays in fact the leading role. 

In the standard uncoupled scenario dark matter perturbations behave as in case [T] above 
(with j4 > and i3 < 0), while dark energy ones provide an example of behavior as in caseEJ 
For coupled models, we concentrate on the case in which the dark-coupling terms dominate over 
the usual one in order to put forward the presence of non adiabatic instabilities (see RefP for 
more details). This strong coupling regime can be characterized by 
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which also guarantees that the interaction among the two dark sectors drives the non-adiabatic 
contribution to the dark energy pressure wave, see Eq. ([5]). At large scales and early times, it 
can be worked out that the main contributions of 5de and 5'^^ coefficients to the second order 
differential equation reduce to: 
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The sign of the coefficient B^ of 6'^^ in this expression is crucial for the analysis of instabilities. 
Assuming c^^^ > 0, it reduces to the sign of the doom factor d defined in Eq. Q. 

Similar second order differential equations for 5de were obtained in Ref. ^ for particular 
expressions of the dark coupling Q and an analytical form of their solutions where derived in 



"Obviously, for \B\ » \A\ a negative B would prevent the onset of growth for any sign of A. 



order to determine when 5de blows up. In particular, the results of Ref.'^ confirm those of Ref.^ 
for positive Q oc pdm and 1 + > 0. In comparison, our approach gives rise to general conditions 
to avoid instabilities as a function of the sign of the doom factor independently of the exact form 
of the dark coupling Q. Indeed, as previously argued, a positive d acts as an antidamping source 
in the growth Eq. ([8|). Whenever d > 1, the overall sign of the A^, coefficient of Sde-, resulting 
from the last two terms in Eq. ([8]), is also positive and it triggers an exponential runaway growth 
of the dark energy perturbations. Large scale instabilities arise then and the universe appears 
to be non viable. 



3 A viable model: Q oc pde 

We have developed a method to determine if a coupled model satisfying Eq. ^ and @ suffer 
or not from non-adiabatic instabilities. We can now easily verify that for 

Q = inpde, (9) 

we have a rather simple and viable model for specific combination oi 1 + w and of the dimen- 
sionless constant coupling ^. In this model, the doom factor of Eq. ([6]) is given by: 

^ (10) 
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Its sign defines the (un)stable regimes. When d < 0, that is, for ^ < and 1 + > (or ^ > 
and \ + w < 0), no instabilities are expected. On the contrary, when ^ and 1 + w have the same 
sign, instabilities will develop at early times whenever d > 1 . 

In Ref. ^ we confirmed these results numerically. Also notice that they are in agreement 
with those of Ref.^, which first pointed out that coupled models with Q = ^Tlpde can be stable 
for 1 + w < 0. They restricted though their stability analysis to the ^ > case. 

In the following, we confront the model satisfying to Eq. ([9]) to cosmological data restricting 
ourselves to negative couplings and w > —1. This guarantees that instability problems in the 
dark energy perturbation equations are avoided for all valued of ^. 



4 Cosmological Constraints from data for Q = ^Tipde 

We explored the current constraints on the dark energy-dark matter coupling ^ using the publicly 
available package cosmomc^. The latter was modified in order to include the coupling among 
the dark matter and dark energy components. More details on the cosmological model and on 
the priors adopted can be found in Ref.'^. The datasets which were taken into account in the 
analysis are: 

1. WMAP 5-year datal^IH! 

2. prior on the Hubble parameter of 72 it 8 km/s/Mpc from the Hubble key project (HST)El 

3. Super Novae (SN) d&ii^ 

4. H{z) data at < 2 < 1.8 from galaxy ages^H^ 

5. large scale structure data (LSS data) from the Sloan Digital Sky Survey^ 

''For Q = ^Tipde, the dark energy density is always positive, all along the cosmic evolution and since its initial 
moment. To ensure that the same happens with the dark matter density, all values of «; < are acceptable for 
^ < 0, while for positive ^ it is required that ^< — w. 



Figure 1: Scenario with Q oc pde- Left (right) panel: la and 2a marginalized contours in the ^-fldmh^ (^~fv) 
plane. The largest, green contours show the current constraints from WMAP (5 year data), HST, SN and H{z) 
data. The smallest, red contours show the current constraints from WMAP (5 year data), HST, SN, H{z) and 

LSS data. 

The data analysis was carried out into two runs, the first run includes the datasets from 1 to 4 
while in the second run the fifth dataset is added. 

Figure [1] (left panel) illustrates the 1 and 2a marginalized contours in the ^-VLdmh^ plane 
where is today's ratio between dark matter energy density and critical energy density. The 
results from the two runs described above are shown. Notice that a huge degeneracy is present, 
being ^ and ^dmh'^ positively correlated. The shape of the contours can be easily understood 
following our discussion in section [TJ In a universe with a negative dark coupling ^, dark matter 
redshift faster. As a consequence, the matter content in the past is higher than in the standard 
uncoupled scenario for a fixed dark matter density today. The amount of intrinsic dark matter 
(which is directly proportional to Qdmh'^) needed to reproduce the LSS data should decrease as 
the dark coupling becomes more and more negative. We also see that LSS data in the second 
run (red contours in Fig. [1]) give rise to the most stringent constraint on the coupling ^. 

The right panel of Fig. [1] shows the correlation among the fraction of matter energy-density 
in the form of massive neutrinos and the dark coupling ^. The relation between the neutrino 
fraction used here and the neutrino mass for Ni, degenerate neutrinos reads 

_ ^li.h'^ _ Y^fny 1 _ N^my I 

~ Qdmh^ ~ 93.2eV ' Qdmh^ ~ 93.2eV ' ^dmh^ ' ^ ' 

Neutrinos can indeed play a relevant role in large scale structure formation and leave key signa- 
tures in several cosmological data sets. Specially, non-relativistic neutrinos in the recent Universe 
suppress the growth of matter density fluctuations and galaxy clustering. This effect can be com- 
pensated by the existence of a coupling between the coupled sectors, given that in this model 
negative couplings enhance the growth of matter density perturbations. 

5 Conclusion 



In this talk, we discuss the origin of non-adiabatic instabilities in coupled models satisfying to 
Eqs. ([3]) and dH). We show that the sign of the doom factor d 



which is a function of the dark couphng Q characterizes the (un)stable regime. In particular, 
when d is positive and sizeable, d > 1, the dark-coupling dependent terms may dominate the 
evolution of dark energy perturbations, which will then enter a runaway, unstable, exponential 
growth regime. 

In a class of viable model in which Q = ^Ji-Pde we have then studied the constraints from 
cosmological data on the dimensionless coupling ^. This analysis was carried out in the ^ < 
and positive (1 + w) region of the parameter space which offers the best agreement with data on 
large scale structure formation. Both w and ^ are not very constrained from data, and it can be 
shown^that substantial values for both parameters, near -0.5, are easily allowed. Furthermore, 
^ turns out to be positively correlated with Qdm.h'^ and larger neutrino fraction is allowed for 
negative ^. 

Acknowledgments 

The work reported has been done in collaboration with B. Gavela, D. Hernandez, O. Mena and 
S. Rigolin and the author was partially supported by CICYT through the project FPA2006- 
05423, by CAM through the project HEPHACOS, P-ESP-00346, by the PAU (Physics of the 
accelerating universe) Consolider Ingenio 2010, by the F.N.R.S. and the I.I.S.N.. 

References 

1. Jussi Valiviita, Elisabetta Majerotto, and Roy Maartens. Instability in interacting dark 
energy and dark matter fluids. JCAP, 0807:020, 2008. 

2. M. B. Gavela, D. Hernandez, L. Lopez Honorez, O. Mena, and S. Rigolin. Dark coupling. 
2009. 

3. Antony Lewis, Anthony Challinor, and Anthony Lasenby. Efficient computation of CMB 
anisotropies in closed FRW models. Astrophys. J., 538:473-476, 2000. 

4. Hideo Kodama and Misao Sasaki. Cosmological Perturbation Theory. Prog. Theor. 
Phys. SuppL, 78:1-166, 1984. 

5. Jian-Hua He, Bin Wang, and Elcio Abdalla. Stability of the curvature perturbation in 
dark sectors' mutual interacting models. Phys. Lett., B671: 139-145, 2009. 

6. Antony Lewis and Sarah Bridle. Cosmological parameters from CMB and other data: a 
Monte- Carlo approach. Phys. Rev., D66:103511, 2002. 

7. J. Dunkley et al. Five- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observa- 
tions: Likelihoods and Parameters from the WMAP data. Astrophys. J. SuppL, 180:306- 
329, 2009. 

8. E. Komatsu et al. Five- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observa- 
tions: Cosmological Interpretation. Astrophys. J. SuppL, 180:330-376, 2009. 

9. W. L. Freedman et al. Final Results from the Hubble Space Telescope Key Project to 
Measure the Hubble Constant. Astrophys. J., 553:47-72, 2001. 

10. M. Kowalski et al. Improved Cosmological Constraints from New, Old and Combined 
Supernova Datasets. Astrophys. J., 686:749-778, 2008. 

11. Joan Simon, Licia Verde, and Raul Jimenez. Constraints on the redshift dependence of 
the dark energy potential. Phys. Rev., D71:123001, 2005. 

12. Max Tegmark et al. Cosmological Constraints from the SDSS Luminous Red Galaxies. 
Phys. Rev., D74: 123507, 2006. 



